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Abstract 

We determine the value of the quark condensate from quenched QCD simulations 
on the lattice in two ways: (i) by using the Gell-Mann-Oakes-Renner (GMOR) for- 
mula; (ii) by comparing the OPE prediction for the Goldstone pole contribution to the 
pseudoscalar vertex, at moderately large momenta. In the MS scheme at fj, = 2 GeV, 
from the GMOR formula we obtain (qq) = - (273 ± 19 MeV) 3 . We show that the 
value extracted from the pseudoscalar vertex, (qq) = — (312 ± 24 MeV) 3 , although 
larger, is consistent with the result obtained from the first (standard) method. 



1 Motivation 



The value of the quark condensate was, and still is, a subject of some controversies. It 
has been experimentally established that in the theory with the spontaneous symmetry 
breaking pattern SU(2)l <S> SU(2)r — > SU{2)y, the quark condensate is indeed the order 
parameter [1]. The extension to the three flavour case still needs to be clarified (for a 
recent critical discussion see ref. [2]). Lattice QCD provides, in principle, the method 
for determining the value of the quark condensate and for studying its dependence on the 
number of dynamical quark flavours. Up to now, the determination of the quark condensate 
on the lattice was limited to the quenched QCD (i.e. with n F = 0). 2 Before tackling the 
theory with n-p = 2 and np = 3 flavors, one would like to learn as much as possible from 
the quenched theory. For example, one would like to understand if the values of the chiral 
condensate obtained by using different methods are consistent among themselves. 

^nitc mixte de Recherche du CNRS - UMR 8627. 

2 See ref. [3] for recent results and the exhaustive list of references. 



The standard method relies on the use of the GMOR formula, i.e. on the same set of 
the background gauge field configurations one computes both the quark masses (m q ) and 
the corresponding pseudoscalar meson masses (mp), and from the slope 

mp = 2B m q , (1) 

one gets an estimate of the quark condensate as 

B Q = -j 2 {qq), (2) 

where / is the pion decay constant in the chiral limit. The renormalization scale and 
scheme dependence of the chiral condensate is just the inverse of the one for the quark 
mass which was discussed in great detail in ref . [4] . 

An alternative way for extracting the value of the quark condensate is from the study 
of the amputated pseudoscalar vertex function, where the cpi^q operator is inserted at 
momentum zero. At moderately large p 2 (p being the momentum flowing through the legs 
of the vertex function), one can compare the shape of this function with the corresponding 
expression derived by means of the operator product expansion (OPE), in which the quark 
condensate appears in the coefficient of the leading power correction. The lattice estimate 
based on this strategy, which is the purpose of this letter, has not been presented so far. 
We show that its value in the continuum limit is fully consistent with the standard value, 
obtained by using the GMOR formula, whose value we updated here as well. 



2 Pseudoscalar vertex 

In this section we discuss the relation between the pseudoscalar vertex and the quark 
propagator and study the dependence of these functions on the chiral quark condensate, 
which enters their OPE as a leading power correction. 

The starting point is to define the quark propagator and the Green function of the 
pseudoscalar density with zero momentum insertion, 

S(p) = Jdx e-*"(q(x)q(0)) > G P (p) = J dx dy e"^-%(:r) g(0) 7 5?(0) ?(y)> ■ (3) 
The amputated vertex function, 

A P (p) = S- 1 (p)G P (p)S' 1 (p), (4) 
is then conveniently projected onto its tree level value 

Tp(p) = ±Tr[ l5 A P (p)} , (5) 

where the trace goes over Dirac and color indices so that the factor 1/12 simply provides 
the normalization to unity. 
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If we write the bare (lattice regularized) inverse quark propagator as 



(6) 



then the basic RI/MOM renormalization condition for the quark propagator in the chiral 
limit can be written as [6] 3 
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where Z q {p) is the quark field renormalization (S(p,p) = Z q (p,)S(p)). 

By studying the quark propagator at large momenta, one can get an estimate of the 
quark mass value, in the RI/MOM scheme, as 
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This estimate has been already discussed in ref. [8]. At lower momenta, however, this 
definition of the quark mass suffers from the presence of the long distance contributions 
due to the coupling to the Goldstone bosons. 

The effect of the Goldstone boson is more clearly seen by considering the quark Ward 
identity which relates the inverse quark propagator to the amputated pseudoscalar Green 
function, 



l5 S-\p 2 ) + S~V)75 = 2Z AP A P (p 2 



(9) 



where Zap is the quark mass obtained from the hadronic axial Ward identity on the 
lattice. 4 After multiplying eq. (9) by 75 and by taking the trace of both sides, we have 



S 2 (p 2 ) = Z AP T P (p 2 ) . 



(10) 



For light quark masses and moderately large momenta, the vertex function Tp(p 2 ) is af- 
fected by the long distance effects which are due to the presence of the Goldstone boson [5] , 
which by means of the LSZ reduction formula generates the term proportional to the Gold- 



stone boson propagator l/(g z 



mt 



Since the operator is inserted at zero momentum, 



q = 0, the vertex function in the chiral limit developes a pole oc l/m 2 ~ 1/p. To account 
for that effect, we expand the vertex function in powers of the quark mass, 



r P (^ P ) = r- btr V) + 



Zap 



+ C(p 2 )p, 



(11) 



3 In practice, we are away from the chiral limit, but the renormalization condition applies equally well 
for m 2 q /p 2 < 1 [7]. 

4 Rccall that 2p = a o (^A o (x)P(0))/(^P(a;)P(0)), with P = q l5 q, A = qj j 5 q, and Z A = Z A (g 2 ) 

X X 

is the (known) axial current renormalization constant. 
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where the first term is the subtracted pseudoscalar vertex, from which the hadronic (Gold- 
stone boson) contribution oc l/m 2 p oc 1/p is subtracted away. The third term is the linear 
quark mass correction while the higher order terms in the expansion, as well as the logarith- 
mic quark mass dependence, are neglected since we deal with light quark masses varying 
in a short interval. 

The renormalization constant of the pseudoscalar density, Z^ MOM (p), is defined in 
terms of the subtracted Green function of eq. (11) through the RI/MOM renormalization 
condition 



-psubtr. 
1 P 



1. 



(12) 



As we already discussed in ref. [4], the value of Zp l ^ MOM (p) obtained from eq. (12) is 
completely consistent with the one obtained by applying the method of ref. [9], which 
allows one to circumvent the second term on the r.h.s. of eq. (11). 

After inserting eq. (11) in (10), multiplying both sides by Z~ 1 ( / u), and accounting for 
the renormalization condition (12), we have 



Mp 2 



Si(p 2 ) 



Za P , B(p 2 



(13) 



RI/MOM 
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where contributions quadratic in the quark mass have been neglected. The first term on 
the right-hand side is the usual short distance quark mass, renormalized in the RI/MOM 
scheme, derived from the axial Ward identity. Equation 13 differs from eq. (8) for the pres- 
ence of the second term on the r.h.s., which represents the power suppressed contribution 
coming from the Goldstone boson. It has been shown long ago that, at the leading order 
in the OPE, this term has the form [10] 



B( P 2 ) 



Mp 2 ) 



OPE 



ctf,,)« + 0(l//) 



(14) 



From this relation we will derive our first estimate of the quark condensate. 

The Wilson coefficient, c(p 2 , p), has been computed at the next-to-leading order (NLO) 
in QCD perturbation theory [11]. In the MS scheme, by choosing the Landau gauge (in 
which the lattice calculations are most easily made), and after setting p 2 = p 2 , one has 5 



c M V) = ~a 8 (p) 



An 



1 + 



99 

T 



10 

n-p 

9 



a sip) 
4ir 



(16) 



We notice that the radiative corrections are large so that at moderately large p 2 they must 
be included in the analysis when extracting the value of the condensate from the lattice 



3 For completeness, we recall the expression for the 2-loop running coupling 
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data. Besides, the inclusion of the radiative corrections is also necessary for specifying the 
renormalization scheme (the leading order anomalous dimension of the quark condensate 
is universal for all renormalization schemes). To eliminate the scale dependence of the 
condensate one defines the renormalization group invariant (RGI) quark condensate, which 
at NLO in perturbation theory is related to the MS one through 



(qq) MS (p) = (<* s (p))-*X 
4 
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3 Lattice data and extraction of the quark condensate 

We work with the 0(a) improved Wilson quark action and use the data-sets consisting of 
(9(1000) independent gauge field configurations, obtained at four different lattice spacings, 
corresponding to (5 — 6.0, 6.2, 6. 4, and 6.45. More complete information about the data- 
sets, as well as the improvement coefficients with the appropriate list of references can 
be found in refs. [4,12]. Since we work at four different lattice spacings, we are able 
to extrapolate to the continuum limit. To eliminate the lattice spacing from the results 
obtained at each lattice coupling, we use the ratio a/r computed in ref. [13], 

(a/ro)p = {0.1863 6 . , 0.1354 6 . 2 , 0.1027 6 . 4 , 0.0962 6 . 45 } , (19) 

so that all our results will be expressed in units of the scale r Q . To convert into physical 
units we will use r = 0.530(25) fm, which corresponds to a^ 60 = 2.0(1) GeV. We will 
also need the quenched value of Aqcd, for which we take r A^g-° = 0.602(48) [14]. 6 

3.1 (qq) from the pseudoscalar vertex 

In order to determine the chiral condensate from the long distance behavior of the pseu- 
doscalar vertex, we first need to extract the function B(p 2 ). That is made by using 10 
different vertex functions, 4 of which are computed with the external legs degenerate in 
the quark mass, and 6 nondegenerate. With these 10 points, for each p 2 , we fit the data 
to the form (11), which we rewrite as 

r P (p 2 , Pi , Pj ) = r- btr V) + z^Tp~) + c{p2){pi + Pi) ■ (20) 
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0.1 0.2 0.3 0.4 

r o(Pi+P 2 ) 

Figure 1: Illustration of the fit to the form (20) at f3 = 6.2, from which we could extract the function 
B(p 2 ), needed for the determination of the chiral condensate. 



The illustration of this fit is provided in fig. 1 for four values of p 2 . We see that the presence 
of the Goldstone pole is indeed pronounced at moderately large values of p 2 . 

Once we identify the Goldstone contribution to the psudoscalar vertex, we perform a 
number of fits to the form 

!^ =< *» w <*^ + 7 + , +Vi (21) 

Li(jr) p z 

where the first term on the r.h.s. is the one that we are interested in (the coefficient c is 
defined in eq. (18)), the second term is the subleading power correction, while the last two 
terms take into account possible contributions of lattice artifacts. To make use of the OPE 
formula we should aim at working at sufficiently large p 2 so that higher powers in 1/p 2 are 
sufficiently suppressed. To do so, we fit the lattice data starting from p cut « 2 GeV, which 
corresponds to (V p cu t) 2 ~ 25, for which the radiative correction term in c RGI (p) is below 
35%. If we set 7 = 5 = A = in (21), then for all our lattices we have xVd.o.f. > 2. 
Therefore one has to let free at least one more parameter. The result of such a fit with 
7 7^ is presented in table 1 and denoted as fit I (see also fig. 2 for illustration). At fixed 
lattice spacing, however, the lattice artifacts may be significant. To examine their impact 
on the value of the quark condensate, we repeat the fits by including either the term with 
p 2 (oc A) or the constant one (oc 5). Both sets of results are reported in table 1, labelled as 
fit II and fit III, respectively. Finally, if we set p cut ^ 3 GeV, the fit with 7 = 5 = A = 
gives a satisfactory x 2 /d.o.f. The corresponding results are denoted as fit IV in table 1. 
We also tried to fit with all the parameters in eq. (21) free (fit V in table 1). 

6 In physical units, A^L =0 = 0.225(20) MeV. 
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Figure 2: Fit of the lattice data to the form (21) for all four lattice spacings considered in this letter. 
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Fit 




-r [(gg) RGI ] 1/3 

U L \ ± ± 1 J 


U ' 


-r 5 x 10 3 


-r„ l \ x 10 6 

U 




6.0 


0.71 ± 0.01 


23 ± 3 






1 


6.2 


0.72 ±0.01 


18 ± 1 








6.4 


0.70 ±0.02 


22 ±3 








6.45 


0.75 ±0.04 


31 ±9 








6.0 


0.82 ± 0.03 


11 ± 5 




29 ± 6 


II 


6.2 


0.74 ±0.02 


17±2 




3.4 ±3.7 




6.4 


0.73 ±0.02 


20 ±3 




4±2 




6.45 


0.73 ±0.04 


34 ± 12 




-4 ±4 




6.0 


0.90 ± 0.04 


5 ± 6 


5 ± 1 




III 


6.2 


0.75 ±0.03 


16 ±3 


1 ± 1 






6.4 


0.75 ±0.02 


18 ±3 


1 ± 1 






6.45 


0.71 ±0.05 


35 ± 14 


-1 ± 1 






6.0 


0.80 ±0.02 








IV 


6.2 


0.80 ±0.01 










6.4 


0.79 ±0.02 










6.45 


0.85 ±0.06 










6.0 


1.14 ±0.10 


-22 ± 15 


27± 11 


112 ± 73 


V 


6.2 


0.90 ±0.07 


5±6 


9±4 


53 ±27 




6.4 


0.85 ±0.07 


10 ±4 


6±4 


24 ± 18 




6.45 


0.87 ±0.06 


24 ± 11 


5±4 


27 ± 19 



Table 1: Details of the fit of the lattice data to the form (21). Various fit forms (labelled as I, 11, III, IV 
and V) are discussed in the text. 

The remaining step towards the determination of the quark condensate is the extrap- 
olation to the continuum limit. Since our action and the renormalization constants are 
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(9(a)-improved, we may attempt extrapolating quadratically in the lattice spacing, i.e. 

roKm)^ 3 = 00 + ^/^)1. (22) 

Our result for Cq, the chiral condensate in the continuum limit, for all four fit forms 
discussed above, are 

r [(^> RGI ] C 1 / 1 = {-0.721(23)i,-0.681(28) n ,-0.672(34) m ,-0.792(24)iv,,-0.54(15)v} -(23) 

One may argue that terms of 0(a) may still be present, since the function B(p 2 ) is obtained 
from the off-shell vertex functions Tp(p 2 , p^, pj), for which the on-shell 0(a) improvement 
does not apply. However, the function B(p 2 ) refers to the chiral limit, and terms in T P (p 2 , p) 
proportional to the quark mass are already taken care of in the fit to the form (11). In 
addition, it has been shown in appendix of ref. [4] that, for these correlation functions, the 
0(a) contribution of operators which are either non gauge-invariant or vanish on-shell by 
the equation of motion vanish in the chiral limit. Therefore, while the 0(d) effects may 
affect the functions r subtr -(p 2 ) and C(p 2 ) in eq. (20) when away from the chiral limit, the 
function B(p 2 ) is polluted by the artefacts 0(a 2 ) and higher. This brings us back to the 
continuum extrapolation form (22). 

What do we learn from the results (23) in the continuum limit? As it can be seen 
from table 1, the corrections oc l/p A are large and positive for every (3. Their neglect in 
the fit IV then expectedly lead to an overestimate of the value for the chiral condensate, 
as confirmed by the last number in eq. (23). Fits II and III give quite consistent values 
for the condensate (in the a — > limit). In other words, the quark condensate in the 
continuum limit is very weakly sensitive to the form of the artifacts that we include in our 
fits (constant or oc p 2 ). The tendency of the artifacts, upon their inclusion in the fit, is to 
lower the value of the condensate. The same tendency is observed also in the fit form V, 
although with larger error bars. 

As our final value we will quote the result of the fit I. The difference between the central 
value of that and the fits obtained by including the artifacts (II and III) is included in the 
systematic uncertainty. The result of the fit V has larger errors and is consistent with the 
results obtained by other fits. As we already pointed out, the radiative corrections are large 
and we take them into account when fitting the lattice data to eq. (21). To account for 
the systematics induced by the omission of higher order corrections in a sip), we will add 
±10% of uncertainty (which represents the square of the 30% effect of the known radiative 
corrections at p = 3 GeV). Finally we have 7 

(gg) RGI = - (269 ± 9tT 8 ± 12 MeV) 3 ± 10% 
<£> (gg) RGI = - (260 ± 9 ± 9 ± 12 MeV) 3 ± 10% 

(gg) M ^(2GeV) = -(312±ll±ll±15±10MeV) 3 , (24) 

where the errors are respectively: statistical, systematics due to the continuum extrap- 
olation, to the uncertainty in r and to the uncertainty due to N 2 LO corrections in the 

7 We remind the reader that r = 0.530(25) fm, is equivalent to r = 2.68(13) GeV -1 . 
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Wilson coefficient c l (p) [see eq. (18)]. Notice that in the second line we symmetrised the 
systematic error bars. 

Finally, we repeated the entire exercise by using the alternative quark mass definition, 
namely the one derived from the vector Ward identity, m q = \ (l/n q — 1/ K crit ), instead of 
the quark mass Zap, used above. The value we obtain in this way is barely distinguishable 
from the one we quoted in eq. (24). 8 

3.2 (qq) from the GMOR formula 



p 


rof 


-r [(gg) RI / MOM (/ia= 1)] 1/3 


-r [(gg) RGI ] 1/3 


6.0 


0.360(7) 


0.709(10) 


0.611(9) 


6.2 


0.365(10) 


0.727(14) 


0.612(12) 


6.4 


0.358(11) 


0.730(16) 


0.605(13) 


6.45 


0.365(39) 


0.743(50) 


0.613(41) 


oo 


0.362(13) 




0.601(25) 



Table 2: Results of the pseudoscalar decay constant in the chiral limit (/) and the chiral condensate, 
obtained by means of the GMOR formula (see eqs. (1) and (2) for all four lattice spacings. We also present 
the results of the linear extrapolation in a 2 to the continuum limit (a — ► 0) . 



We now repeat the standard exercise of extracting the value of the quark condensate 
by employing the GMOR formula. The values of the pseudoscalar meson and the quark 
masses are all listed in table 2 of ref. [12]. In table 2 of the present letter, we give the 
results obtained by using eqs. (1) and (2), where we use for the quark mass the one 
defined via the axial Ward identity (p). The needed renormalization constants, Za and 
Zp^ MOM (l/a), are given in ref. [4]. For completeness, we also present the values of the 
(improved) pseudoscalar meson decay constant in the chiral limit, /, which is obtained by 
linearly extrapolating in the quark masses (fp — f + const, x p). To convert the quark 
condensate from the RI/MOM scheme to the RGI form, we use the anomalous dimension 
known up to 4-loops [15]. These latter results are then extrapolated to the continuum limit 
linearly in a 2 [see eq. (22)]. In physical units, our results read 

(gg) RGI = - (224 ± 9 ± 10 MeV) 3 (qq) m (2 GeV) = - (273 ± 11 ± 15 MeV) 3 . (25) 

We checked that this value is completely consistent with the alternative definition of the 
quark mass, namely with m q = | — 1/ n crit ), and with Z^ MOM (l/a) also given in 

ref. [4]. 9 Finally we also note that the above result agrees very well with the QCD sum 
rule estimate of ref. [16], where (qq) m (2 GeV) = - (267 ± 16 MeV) 3 has been quoted. 

8 More specifically, with m q instead of p, we get (qq) WS (2 GeV) = - (313 ± 11 ± 13 ± 15 ± 10 MeV) 3 . 
9 Whcn using the quark mass defined via the vector Ward identity instead of the axial one, we get 
{qq) m (2 GeV) = - (268 ± 13 ± 15 MeV) 3 . 
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4 Summary and conclusion 



We now briefly summarize our findings and comment on some results reported in the 
literature. 

(1) We update the value of the chiral condensate obtained from quenched QCD on the 
lattice, by using the GMOR formula. After combining the errors given in eq. (25) in the 
quadrature, we have 

(99>gmor(2 GeV) = - (273 ± 19 MeV) 3 . (26) 

This result is obtained from simulations performed at four lattice spacings, by employing 
non-perturbative renormalization and C(a)-improvement, followed by an extrapolation to 
the continuum limit. 

(2) We compute the quark condensate by using an alternative strategy, namely by study- 
ing the long distance (Goldstone) part of the pseudoscalar vertex function. In terms of 
the OPE, the chiral condensate appears in the coefficient of the leading power correc- 
tion in 1/p 2 . From the calculations at four lattice spacings and after extrapolating to the 
continuum limit we obtain 

(99>ote( 2 GeV) = - (312 ± 24 MeV) 3 . (27) 

(3) From the above results, it seems that the two completely different strategies lead to 
quite a consistent value of the quark condensate. To better appreciate this point we rewrite 
the RGI results in units of r , i.e., 

-r [(qq) RGl } 1/3 = {0.701(23)(20)(23) O pe, 0.601(25) GMO r} • (28) 

The method based on using OPE is less reliable since radiative and further power correc- 
tions are large. Even if we combine the errors in quadrature the agreement would be at the 
2cr-level, which is far from what has been claimed in ref. [5], where the OPE and GMOR 
results were argued to differ by a factor of 3. 

Before closing this letter, we should explain why our conclusion is qualitatively different 
from the one reported in ref. [5]. The first difference is that in eq. (11), besides the 
Goldstone term (oc 1/p) we also allow for the presence of the term linear in quark mass. 
Such a term could not be studied in ref. [5] since only three quark masses were considered. 
The net effect of this modification is that the function B(p 2 ) becomes smaller. Secondly, 
in the OPE, we allow for the presence of the term oc l/p A , and we find that, for moderately 
large momenta, this (subleading) power correction is not negligible, while in ref. [5] such 
a term was not found. Finally, we also accounted for the terms that are due to the lattice 
artifacts [see eq. (22)], which further reduce the value of the condensate in the continuum 
limit. Such effects were not studied in ref. [5] (where they only considered the data produced 
at f3 — 6.0). Notice also that the reference value of the chiral condensate considered in 
ref. [5], was 20% smaller than the one we obtain here by using the GMOR formula. 
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